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Abstract. Supernovae of Type Iln (narrow line) appear to be explosions that had strong mass 
f^ ] loss before the event, so that the optical luminosity is powered by the circumstellar interaction. 
T-H . If the mass loss region has an optical depth > c/vs , where Vs is the shock velocity, the shock 
^^ ■ breakout occurs in the mass loss region and a significant fraction of the explosion energy can be 
Cn , radiated. The emission from the superluminous SN 2006gy and the normal luminosity SN 2011ht 
$^ ■ can plausibly be attributed to shock breakout in a wind, with SN 2011ht being a low energy 
C^' event. Superluminous supernovae of Type I may derive their luminosity from interaction with a 
•^r , mass loss region of limited extent. However, the distinctive temperature increase to maximum 

luminosity has not been clearly observed in Type I events. Suggested mechanisms for the strong 

mass loss include pulsational pair instability, gravity-waves generated by instabilities in late 

burning phases, and binary effects. 
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1. Introduction 



^■ 

Oh! 

JL ■ Once the radiation from the supernova shock wave emerges from the steUar surface, 

^ I the interaction with the surrounding medium commences. In the case of core collapse 
[/3 ■ supernovae, the initial interaction is with mass loss from the progenitor star. At low 
C^ [ circumstellar densities, the emission from the interaction is a small part of the super- 
nova power and typically involves X-ray and radio radiation. This is the case for Type 
7-H IIP supernovae, where the inferred circumstellar density is roughly consistent with that 

K«*" ■ expected from the red supergiant progenitors of the these events (Chevalier et al. 2006). 
^^ ' However, there are supernovae, of Type Iln, that are thought to have most of their 
peak optical power come from interaction with a circumstellar medium (e.g., Chugai & 
Danziger 1994). Here I discuss cases where the circumstellar interaction is an important 
part of the optical emission. The shock breakout phase is treated in Section 2 and the 
initial viscous shock propagation in Section 3. Speculations on the formation of the dense 
circumstellar medium are discussed in Section 4 and the final comments in Section 5. 



.!^ 2. Shock Breakout 

K% A useful way of considering the enhanced supernova luminosity that can result from 

S [ very dense circumstellar interaction is in terms shock breakout in the dense gas (Cheva- 

" " " lier & Irwin 2011, Balberg & Loeb 2011). While a shock front is inside the supernova 
progenitor, the shock front is mediated by radiation for standard supernova conditions. 
The optical depth across the shock transition for a radiation dominated shock is ~ c/Vs, 
where Vg is the shock velocity and c is the speed of light. This result is determined by the 
condition that the hydrodynamic time for the shock wave to cross the shock thickness is 
approximately equal to the diffusion time for the radiation. When the shock front reaches 
a position that is at an optical depth of ^ cjvg from the surface, radiation can diffuse 
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out to optically thin layers and shock breakout occurs. For low circumstcUar densities, 
the timescale for the breakout event is typically quite brief and depends on the radius 
of the progenitor star, R. For a red supergiant progenitor, the most extended normal 
progenitor, the supernova shock takes about a day to traverse the star and the breakout 
phase lasts about half an hour, determined by the light travel time across the star, R/c. 

A low density, optically thin circumstellar medium is not expected to have much in- 
fluence on the optical light near maximum, although it can be very important for radio 
and X-ray emission. At higher circumstellar densities, the medium becomes optically 
thick and the breakout radiation must diffuse out through the circumstellar medium. 
The result is that the radiated energy in the breakout is not changed, but the radiation 
is emitted over a longer period of time, lowering the luminosity. 

When the circumstellar optical depth, r^, becomes of order c/vs or more, the radiation 
dominated shock can propagate into the circumstellar region. In this case, substantial su- 
pernova energy can be present as internal energy at a larger radius than would otherwise 
be the case, so that the effects of adiabatic expansion are lessened and the luminosity 
rises. Chevalier & Irwin (2011) suggested that this model can describe the initial rise of 
lO's of days of superluminous supernovae like SN 2006gy. When the circumstellar inter- 
action occurs at a large radius where diffusion of radiation can occur rapidly and the 
circumstellar medium is sufficiently massive to thermalize the supernova kinetic energy, 
a large fraction of the supernova energy can be radiated. The physical parameters, a cir- 
cumstellar medium of ~ 10 Mq extending out to ~ 10^^ cm, are similar to those in the 
diffusion model of Smith & McCray (2007). Chevalier & Irwin (2011) used an analytic 
model for the emission, which has limitations; in particular, the initial interaction was 
described by a self-similar solution, but that solution breaks down as the reverse shock 
wave propagates into the supernova ejecta. Ginzburg & Balberg (2012) undertook nu- 
merical simulations that did not have this problem and obtained comparable parameters 
for SN 2006gy. In their simulations, Moriya et al. (2013) found that a density law p oc r~^ 
did not give a good fit to light curve; their favored explosion parameters are mass loss 
as high as 15 M© and energy no more than 4 x 10^^ ergs. 

A general expectation of the shock breakout model is that the initial rise to maximum 
luminosity should be primarily due to heating of the photosphere as the shock wave 
begins to affect the stellar emission. This effect is not clearly observed in the data on SN 
2006gy, although there is some sign of an increasing temperature in the spectra of Smith 
et al. (2010). 

A Type Iln event that apparently does show the increasing temperature in the rise 
to maximum is SN 2011ht. Early Swift observations show a clear rise in temperature 
to maximum (Roming et al. 2012), which is characteristic of shock breakout in the cir- 
cumstellar gas. The rise time to maximum, 40 days, is an indicator of the density of 
the circumstellar gas. The supernova radiation is expected to be efficiently radiated, so 
that the radiated energy is a significant fraction of the supernova energy if the energy is 
thermalized by the interaction.. Interestingly, the light curve is similar in shape to that 
of SN 2006gy, but the absolute magnitude is fainter by 5 magnitudes, with a radiated 
energy ^ 2.5 x 10^^ ergs. The implication is that the explosion energy is smaller by a 
factor up to ~ 10^. The possibility that the event was not a supernova was noted by 
Roming et al. (2012). However, there are supernovae that occur in this low energy range 
and Mauerhan et al. (2013) suggest that the late emission from SN 2011ht is due to a 
small amount of "''^Ni. In the end, it is not possible to definitively answer the question of 
whether SN 2011ht was a supernova. 

For luminous supernovae that are of Type Iln, it is very plausible that the high optical 
luminosity is produced by circumstellar interaction (see next section). There are also 
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luminous supernovae that are not of the "n" type. SN 2008es is a luminous Type II 
(Gezari et al. 2009, Miller et al. 2009) and SN 2010gx is an example of a luminous 
Type Ib/c (Pastorello et al. 2010). The Type Ib/c objects outnumber those of Type II. 
Chevalier & Irwin (2011) suggested that the optical light from these objects is powered by 
circumstellar interaction in which the mass loss region ends before the diffusion radius is 
reached by the shock front. Ginzburg & Balberg (2012) carried out numerical simulations 
of this situation and successfully modeled the optical light curve from SN 2010gx. As 
described above, an expectation in this model is that the temperature should increase in 
the rise to maximum luminosity. In the case of SN 2006oz, there are observations on the 
rise to maximum and there is no evidence for an increasing temperature; the temperature 
remains roughly constant (Leloudas et al. 2012). Another mechanism, such as magnetar 
power (Kasen & Bildsten 2010), may be indicated. 

Another superluminous event of interest is PSl-lOafx (Chornock et al. 2013). In this 
case, the luminosity and temperature at maximum implied a photospheric radius of 
5 X 10^^ cm, while the observed lines indicated an expansion velocity of 11,000 km s^^. 
The implied expansion timescale is 50 days, but the rise time for the light curve is 
~ 10 — 15 days. This evolution is naturally explained in the shock breakout view (Cheva- 
lier & Irwin 2011); the shock front must traverse the optically thick star before shock 
breakout occurs. However, the rise to luminosity maximum did not show the increasing 
temperature that might be expected in the breakout scenario and there is no spectral 
evidence for circumstellar interaction, so the nature of the event is still in doubt. 

Overall, the shock breakout scenario is attractive for the events that show narrow line 
spectra because there is evidence for slowly expanding circumstellar matter ahead of 
the shock front. When these lines are not present, there is no clear evidence for dense 
circumstellar matter. There may be more than one mechanism that gives rise to very 
luminous supernovae. 



3. Viscous Shock Wave 

Once the forward shock wave gets into a regime where the optical depth to the surface 
T < c/vs, a radiation dominated shock can no longer be maintained and there is a 
transition to a viscous shock front. There may initially be radiative acceleration of the 
unshockcd gas, so that a viscous shock does not initially form, but, if the circumstellar 
medium is extended, the formation of a viscous shock in inevitable. The maximum optical 
depth at which it can form is somewhat less than c/vs- 

At low circumstellar densities the emission from the reverse shock dominates that from 
the forward shock because the density is higher and the shock velocity is lower, ^ 1000 
km s^^ or less. However, as the circumstellar density rises, the reverse shock becomes 
radiative first, which leads a lower increase in luminosity with increasing density and to 
a dense shell that can absorb radiation from the reverse shock, thus affecting the X-ray 
luminosity. 

The forward shock, with a velocity of lOOO's of km s""'^ or more, can heat gas to 
a temperature of 100w| KeV, where w| is the shock velocity in units of 10^ km s^^. 
Katz et al. (2011) found that inverse Compton cooling can balance the shock heating 
at a temperature of 60 keV if the optical depth is c/vs (breakout conditions) and the 
shock velocity is 10^ km s^^. This estimate assumed that the electrons are heated only 
by Coulomb collisions with the ions. If there is additional collisionless heating, the gas 
temperature is raised. Chevalier & Irwin (2012) estimated regions of the shock velocity vs 
circumstellar density plane where the cooling time is short compared to the age and where 
inverse Compton cooling is larger than bremsstralung cooling. The regime where the 
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cooling is rapid roughly overlaps the regime where the electron scattering optical depth 
is > 1. If inverse Compton losses are larger than bremsstrahlung, the X-ray emission 
from the hot shocked gas is reduced as a fraction of the shock power. 

The X-rays that are emitted by the hot gas have to propagate through the preshock 
circumstellar medium. One effect of the surrounding medium is Compton recoil (Chevalier 
& Irwin 2012, Svirski et al. 2012), which reduces the energies of the high energy photons 
to a maximum of ^ 511/Tg^ keV, where t^s is the electron scattering optical depth and 
is assumed > 1. It can be seen that Tes > 8 is necessary to have an effect in the energy 
range around 10 keV that is accessible to X-ray telescopes like Chandra. NuSTAR, with 
its sensitivity to 80 keV, has a better chance of detecting this effect. 

Another effect is photoabsorption by the preshock medium (Chevalier & Irwin 2012), 
which has a larger effect at low X-ray energies. An important aspect of photoabsorption 
is that its occurrence requires that the absorbing medium not be completely ionized. For 
the X-ray emission that is expected from the hot gas, there is the possibility of complete 
ionization (Chevalier & Irwin 2012). Complete ionization is expected at higher shock 
velocities ^ 10"* km s~^, but not at lower velocities, < 5000 km s~^. The result also 
depends on clumping and asymmetries in the circumstellar gas. If the gas is not fully 
ionized, the column density corresponding to an electron scattering optical depth Tes is 
3 X 10'^^Tes{Z/ZQ)~^ cm~^, where Z is the supernova metallicity and Zq is the solar 
metallicity. Even for Tes = 1, the column density is orders of magnitude larger than a 
typical interstellar column density and would have a dramatic effect on X-ray emission. 

A supernova that provides some test of these models is the Type Iln (narrow line) SN 
2010J1. The optical spectrum showed narrow, presumably circumstellar, emission lines on 
top of broad wings that could be attributed to electron scattering (Smith et al. 2010). 
The scattering requires an electron scattering optical depth ^ 1. An X-ray spectrum with 
the Chandra observatory at 2 months required a column density ^ 10^** cm^^, (Chandra 
et al. 2012). At an age of 12 months, the column dropped to ~ 3 x 10^^ cm~^, confirming 
that the absorption is connected to circumstellar gas. The observations with Chandra 
showed a hot thermal component with temperature ^ 10 keV, corresponding to a shock 
velocity ^ 3000 km s^^. The emission is presumably from the forward shock front in SN 
2010J1. 

The presence of the large column density in SN 2010J1 shows that it is possible for 
a large fraction of the X-ray emission to be extinguished. The absorption in SN 2010J1 
corresponds to Tes ^ 1; larger values of Tes are expected close to the time of shock 
breakout in a wind. The expected absorption optical depth at 10 keV is ^ \.hTes{Z JZq) 
if the circumstellar region is not fully ionized. Chevalier & Irwin (2012) noted that the 
observed low X-ray luminosity of SN 2006gy, ^^ 10"*" ergs s^^ (Smith et al. 2007; Ofek 
et al. 2007) compared to an optical luminosity of 3 x lO**^ ergs s~^ at a time close to 
maximum light (Smith et al. 2010), is consistent with breakout in a wind. The case of SN 
2011ht was also mentioned above as a likely case of a wind breakout. An X-ray luminosity 
of lO^'' ergs s~^ was initially reported for this object based on Swift observations (Roming 
et al. 2012), but higher spatial observations with Chandra showed that the source is not 
coincident with SN 2011ht (Pooley 2012). Taking a conservative upper limit to the X-ray 
luminosity of 10^^ ergs s^^, the X-ray emission is again a small fraction of the optical 
luminosity, which was 3 x 10**^ ergs s~^. The strong suppression of the X-ray emission in 
SN 2006gy and SN 2011ht is consistent with both of these objects being shock breakouts 
in a wind, although the explosion energy is probably much smaller in SN 2011ht. 

Dwarkadas & Gruszko (2012) make the point that the evolution of X-ray emission 
from supernovae is typically not what is expected for a steady wind, so that estimates 
of mass loss rate can be misleading. However, the preshock density at any time can be 
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estimated from po = L/(2TTriR^Vg), where 77 is the efficiency of conversion of the shock 
power into the luminosity L. If observations are made over a period of time, an estimate 
of the mass required for the luminosity can be made and if the velocity of the mass loss 
can be estimated from narrow lines, the timescale for the emission can be obtained. The 
mass and timescale estimates provide useful constraints on possible mechanisms for the 
mass loss. 

An interesting aspect of the viscous shock front in dense mass loss is that the shock 
wave is expected to be coUisionless, so that diffusive shock acceleration of particles to 
relativistic energies can occur (Katz et al. 2011, Murase et al. 2011). These studies show 
that an unusually nearby dense interaction supernova would have to occur to have a 
chance of detecting high energy gamma-ray emission from the supernova. 

4. Origin of the Dense Mass Loss 

Many of the characteristics of Type Iln supernovae require an optical depth to electron 
scattering ^ 1 in the circumstellar medium, which can be expected to produce wings to 
narrow line profiles and obscure the inner high velocities related to the supernova (if 
present). Typical parameters for Type Iln events involve masses of 0.01 — 10 Mq on 
scales of 10^^ — 10^^ cm from the supernova. As noted above, a circumstellar mass as 
high as 10 — 20 Mq has been estimated for SN 2006gy. These events have been linked 
to luminous blue variables (LBVs) because they are some of the only objects known to 
have such extreme mass loss (Kiewe et al. 2012). In addition, the Type Iln SN 2005gl was 
found to have a probable LBV progenitor in pre-explosion images (Gal- Yam & Leonard 
2009). Although this link has been made, there is still the issue that the reason for the 
extreme mass loss from LBVs is not understood. In addition, stellar evolution models 
do not predict supernova explosions at the time of the LBV phase. The timing of the 
explosion close in time to the mass loss is a remarkable feature of the Type Iln events. 

One explanation for SN 2006gy was pulsational pair instability of massive stars (Woosley 
et al. 2007). In this model, there is no terminal explosion, but the emission results from 
the interaction of pulsationally driven shells. The shells can have a mass of several Mq 
and the radiated energy from their interaction can be 10^" ergs. 

An interesting proposal by Quatacrt & Shiode (2012) is that super-Eddington fusion 
reactions at the end of the life of a massive star can generate gravity waves that deposit 
their energy/momentum in the outer parts of the star, driving strong mass loss. The 
prime burning phases for this are Ne and O burning, which occur in the last year of 
evolution before core collapse. However, only a fraction of massive stars show evidence 
for the late mass loss and, in this hypothesis, it is unclear what determines the fraction 
that show the mass loss. A variant of this hypothesis is that the burning instability 
leads to expansion of the star instead of mass loss, and the mass loss is driven by binary 
interaction (Soker 2013). 

Ofek et al. (2013) found emission from SN 2010mc in the 40 days leading up to the 
explosion and suggested that the observations support the gravity wave driven mass loss 
hypothesis. The reasons are that the circumstellar velocities are relatively high, ^ 2000 
km s^^ and that the brightening was observed 40 days before explosion, when the late 
unstable burning phases are expected. However, Smith et al. (2013) show that the light 
curve of SN 2010mc is very much like that of SN 2009ip during the period June - Sept 
2012, indicating that they involve a similar physical situation. SN 2009ip showed eruptive 
behavior over the period 2009 to 2012, somewhat longer than might be expected for the 
Ne and O burning phases. Also, high velocities observed in spectral lines are possibly 
the result of electron scattering. In any case, Fraser et al. (2013) note that there is some 
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question whether SN 2009ip was in fact a supernova in Sept 2012, a question that can 
also be raised for SN 2010nic. 

Another possibility for the strong mass loss is that it is driven by close binary inter- 
action. Common envelope evolution is expected in close binaries, accompanied by strong 
mass loss. This has the advantage that there can be considerable variety in the events. 
Also, recent observations have shown that massive stars are in close binaries more fre- 
quently than had been thought (Sana et al. 2012). The problem is to have an explosion 
close to the time that the mass loss occurs. Chevalier (2012) speculated that the event 
involves the spiral in of a neutron star due to common envelope evolution. In cases where 
the interaction occurs at a relatively small separation, the neutron star spirals into the 
core of its companion because of the steep density gradient. In cases where there is spiral 
in, but the separation is somewhat greater, the companion star can become a red super- 
giant with a flat density profile, which causes spiral in to stop. The result is a neutron 
star - helium star binary. When the He star evolves, there is again the possibility of mass 
loss accompanied by the spiral in of the neutron star to the core. 

In one view, the spiral in of the neutron star to the core gives rise to a red supergiant 
with a neutron star core, a Thorne-Zytkow star (Thorne & Zytkow 1977). However, 
neutrino losses of the matter near the neutron star may lead to strong accretion onto 
the neutron star, especially when the neutron star is in the core of the companion star 
(Chevalier 1996). The rapid accretion may lead to black hole formation and an explosion, 
as in the model for gamma-ray bursts (GRBs) of Fryer & Woosley (1998). A variant of this 
model is that the rapid accretion is accompanied by strong magnetic field amplification 
of the neutron star, leading to an explosion (Barkov & Komissarov 2011), 

In this scenario for explosions preceded by strong mass loss, the explosion mechanism is 
not the neutrino mechanism, which is generally favored for most core collapse supernovae, 
but is related to rapid rotation and accretion on a central compact object. Spiral in of a 
compact object is an efficient way of achieving rapid rotation in the central region. The 
mechanism is believed to operate for GRBs, which have been associated with energetic (> 
10^^ ergs) supernovae. Some Type Iln supernovae also appear to have had unusually large 
energies; for example, SN 2003ma had an integrated bolometric luminosity of 4 x 10^^ 
ergs over 4.7 years and an estimated explosion energy > 10^^ ergs (Rest et al. 2011). 
Studies of the neutrino mechanism for core collapse supernova explosions indicate that 
the maximum energy that can be attained is ~ 2 x 10^^ ergs (Janka 2012). Both GRBs 
and some Type Iln supernovae may require an explosion mechanism involving rotation 
and magnetic fields in order to produce high energy explosions. 



5. Discussion 

Although circumstellar interaction at low densities can be successfully described in 
terms of spherical models, there are signs that at the higher densities considered here, 
where the interaction contributes to the optical luminosity, the situation is more complex. 
Galactic LBVs show a complex circumstellar structure, binary evolution can lead to 
aspherical strcuture (as apparently occurred in SN 1987A), and Type Iln supernovae 
often show significant polarization in their optical light. In addition to these complexities, 
the Type Iln events appear to include a wide range of stellar masses and explosion 
energies, indicating a variety of evolutionary paths. We are still at the early stages of 
understanding these events. 
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